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Abstract. - The dynamics of sandpile models are mapped to discrete interface equations. 
We study in detail the Bak-Tang-Wiesenfeld model, a stochastic model with random thresh- 
olds, and the Manna model. These sandpile models are, respectively, discretizations of the 
Edwards-Wilkinson equation with columnar, point-like and correlated quenched noise, with 
the constraint that the interface velocity is either zero or one. The constraint, embedded in 
the sandpile rules, gives rise to another noise component. Studies of this term for the Bak- 
Tang-Wiesenfeld model reveal long-range on-site correlations and that with open boundary 
conditions there is no spatial translational invariance. 



Simple cellular automata have been under intensive study during the last decade. The 
aim has been to explain many of the power-laws that can be often seen in systems in nature. 
A particularly famous example is the paradigmatic Bak-Tang-Wiesenfeld (BTW) sandpile 
model jjj , which exhibits both spatial and temporal criticality using the language of statistical 
mechanics. Many variants of sandpile models have been developed, while the BTW model 
still occupies a central position due to its deceptive simplicity and complicated behavior || . 

The central ingredients of a sandpile model are rules for redistributing grains and for 
driving the system by grain addition if it is stable. This gives rise to two time scales, slow 
and fast, and usually to diffusive grain dynamics. Grains are removed from unstable sites and 
distributed to neighbors. Those lost through the open boundaries are compensated by adding 
grains randomly, but on a slow time scale compared to the time that individual avalanches or 
periods of activity take. In the asymptotic self-organized critical (SOC) steady state power- 
laws emerge in the avalanche statistics. The scenario is reminiscent of another problem in 
non-equilibrium statistical mechanics: driven interfaces which couple a random environment 
to an elastic object. At a critical value of the drive force the velocity of the interface becomes 
non-zero, and exactly at this depinning transition critical behavior ensues. The time and 
length scales are renormalized, with associated critical exponents 

In this Letter we demonstrate a connection between interfaces and many sandpile models. 
We construct a mapping of the BTW model, and two stochastic models ("rice-pile" models 
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the Manna- model ||) to discrete models of driven interfaces (9). The interface field is 
constructed such that it counts topplings as function of time. By employing the mapping, one 
obtains a description of the stochastic movement of sandpile grains as an interface propagating 
through a quenched random medium. Such a description of sandpile dynamics allows us to 
address several issues in SOC. First, the universality classes of different models are reflected 
in the noise terms of the interface equation. The quenched noise arises from the shot noise of 
grain deposition, and this is the important feature in the BTW model. The noise can also result 
from the fluctuations in the movement of grains that is caused by randomness in the rules. 
Finally, disorder is created by the fluctuations in the density of active (unstable) sites. A grain 
may sometimes end on a site that would become anyways unstable. Such over-active sites 
correspond to interface locations with a local velocity constrained to unity. This phenomenon 
is most distinctive in the BTW model, as compared to the stochastic threshold/rice-pile 
and Manna models. The rice-pile model is exactly a discrete version of the (random field) 
linear interface model class and may be discussed with well-known scaling and renormalization 
arguments The mapping establishes the upper critical dimension to d u — 4 for all the 

models. As a novel effect, the noise coming from the velocity constraint demonstrates a lack 
of translational invariance in open systems. 

The sandpile models are as follows: each site x of a hyper-cubic lattice of size L d has z(x, t) 
grains. When z(x, t) exceeds a critical threshold z c (x), the site is active and topples. Grains are 
removed from x and given to the nearest neighbors (nn) . This in turn may cause some nn's to 
topple and so on. If there are no active sites in the system, one grain is added to a randomly 
chosen site, z(x,t) — > z(x,t) + 1. The time and number of topplings till the system again 
contains no active sites define an avalanche and its internal lifetime. The system is usually 
open such that grains which topple out of the system are lost (in d = 1: zq = = 0). The 
specific toppling rule distinguishes between the models. The BTW model has z c equal to a 
constant, z c — 2d — 1 and the toppling rule: z(x, t + 1) = z(x, t) — 2d, z(y, t + 1) = z(y, t) + 1, 
where y denotes all the 2c? nn's of the site x. The rice pile model has the same rule but 
with z c (x) randomly chosen after each toppling from a probability distribution. In the Manna 
model, two grains are given to two randomly chosen neighbors with z c — 1. We consider here 
two kinds of ensembles. The SOC one is defined above via the drive and the open boundaries. 
Another possibility is to use periodic boundary conditions, and prepare the system at a certain 
average 'energy' (z(x, 0)) in which case the critical state is reached at a certain energy (z) c flC)|| . 

The mapping of the dynamics of sandpile models ||. begins with the definition of an 
interface or memory field H (x, i). It counts topplings at site x up to time t. The dynamics of 
H defined in this way reads 

where we have rewritten the equation of motion for H in the form of a discrete interface 
equation. The important point is the construction of the local "force" f{x, t) = z(x, t) — z c (x) 
at site x. We can express f(x, t) in terms of the grains added to site x, n™, and removed from 
x, n° ut , as f(x, t) = n™ - n% ut - z c (x), since z(x, t) = n™ - n° ut . 

For the BTW and rice-pile models n™ and n° ut can be derived from the local height field 
H(x,t) and an external force term F(x,t). Here, F counts the number of grains added to 
site x by the external drive up to time t. Thus for SOC drive F(x, i) increases on the slow 
time scale and does not change during avalanches, i.e. it acts as a (columnar) quenched noise. 
The construction of / using F works for any drive (e.g. continuous, uniform). Observing that 
n™' is simply n° ut — 2dH(x,t) and n™ = J2 X H(x nn ,t) + F(x,t), where x nn denotes the 
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2d nn's of x, one arrives at f(x,t) — \7 2 H + F(x,t) — z c (x,H), where V 2 H is the discrete 
Laplacian []ll| , |l2]j . The threshold z c (x,H) depends on H for the rice-pile model and reveals 
that the random thresholds correspond to a quenched force field which is random in x and H, 
and acts on the interface. For the BTW model z c is constant. 

For the Manna model we have to incorporate the randomness in the distribution rule. We 
use a projection technique by writing as the average incoming flux n™ plus a fluctuating 
part 6n l ™ and arrive at 

f(x, t) = l/d V 2 H + F(x, t) - z c (x, H) + t(x, H), (2) 

with the resulting noise term 

r(x,H) = 6n™ = n™ - = n™ - l/d^^x^t). (3) 



Each nn-toppling contributes 1/d to fi l ™ as the sum over the nn's is the expected number of 
grains from neighboring sites due to their topplings. When site x topples, one uses the above 
definition of t(x,H) in terms of the fluxes (n m , n m ) to evaluate it at site x and height H. 
The projection trick, used to construct the quenched noise r, means that the randomness in 
the toppling rule is matched to an equivalent noise field r so that the interface equation for 
H reproduces exactly the behavior of the Manna sandpile. It can also be used in the case 
of other models, where the effect of a toppling is random. Note that t(x,H) is a conserving 
noise since the random toppling rule for the Manna model conserves the number of grains 
(except at open boundaries). 

The step-function, 9(f), in Equation ([!]) forces it so that the interface does not move 
backwards and such that the velocity v = AH /At is either or I ||. This means that 
sandpilcs arc equivalent to cellular automaton models of interface depinning, with this velocity 
constraint. We next map this constraint into an effective noise term, denoted er, in the interface 
equation, that can be used to discuss the possible differences of 'real' depinning models and 
those that arise from sandpile models via the mapping. Consider the toppling example in 
Fig. 0. On the avalanche timescale / < at site x. As function of time, / increases until at 
time t — 1 several neighbors topple resulting in / > 0, so that site x will topple at time t and 
AH/ At = 1. 

The sandpile rules result in an effective force /' = 1 that acts on the interface H at the 
time of toppling: AH/ At = f 9(f) — f'9(f). The relation between / and /' for each toppling 
at x (constant H ) reads thus 

f'(x,t) = f(x,t) + a(x,H) a(x,H) = l + z c (x,H)-z(x,t*) (4) 

where t* is the time at which site x topples such that cr(x, H) by this construction is a quenched 
random variable. It is computed from the difference between / and /' when x topples. Notice 
that /' and / are by definition time-dependent variables, since they change as grains are 
moved, or the Laplacian changes. They also contain a quenched force component, as is seen 
from Eq. (Eh. The easiest way to study the er-noise is to determine it in a simulation using 
Equation (0) . This construction of a is similar to that of the r-noise term for the Manna-model 
and other models for which the projection trick can be used. The trick maps the difference of 
the expected value of the local force at toppling f(x, t*) = z x — z c and the true one to r. Here 
the difference of f (x, t*) = 1 and / maps to a, a quenched variable. Such differences arise in 
the Manna model due to randomness in grain movement, and in the case of the a- noise from 
the effect of the step-function. 
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The point in the noise variables cr(;r, H) and t(x, H) is that given these disorder fields, the 
interface equation exactly reproduces the history or dynamics of a sandpile "run" . Thus we 
can study the interface model as such, and try to infer the properties of the original sandpile 
from its behavior. Also, the cr-noise allows one to study the effects of the peculiar discretization 
(dH/dt = 0(f)) explicitly, since AH/ At = f'O(f') can be interpreted as the discretization of 
the continuum equation dH/dt = /' (5). Combining Eqs. ([[]) and (|J) we write the discretized 
interface equation as 

A TT 

— = uV 2 H + r)(x,H)+F(x,t) + a(x,H). (5) 

Here the diffusion constant v is unity for the BTW and ricepile models and 1/d for the Manna 
model, and the quenched noise rj(x,H) = —z c (x,H) + t(x,H). Equation (|J) is the central 
difference discretization of a continuum diffusion equation with quenched noise, called the 
linear interface model (LIM) or the quenched Edwards- Wilkinson equation j|-|5| . Equation (||) 
contains two main ingredients. First, the Laplacian character of these sandpile models which 
is such that the differences in accumulated topplings map exactly to an elastic force. In 
the interface language once the force increases sufficiently to overcome the pinning force, 
\7 2 H + r](x, H) + F(x, t), the interface moves by one step. For open boundary sandpiles, the 
right interface boundary condition is H = which is to be imposed at "extra sites" (x = 0, 
x — L + 1 for a system of size L in Id). In the SOC steady-state the Laplacian increases, 
because of the roughly parabolic shape for the the toppling profile H(x) (see jll],^0]]). This is 
compensated by the ever-increasing (F(t)}, or the addition of grains by the SOC drive. 

Second, the randomness in the sandpile rules map into noise variables [F, 77], as do the 
details of the dynamics [a]. Equation (||) allows thus to make conclusions about the uni- 
versality classes of models based on the noise terms and their relevance. The F-teim in 
Eq. (||) is columnar: it integrates the deposited grains and is constant during avalanches. The 
T](x, H) term in the rice-pile and Manna models explicitly depends on H . For the rice-pile 
T](x, H) = —z c (x, H), trivially. The associated LIM has point-disorder since z c (x, H) is delta- 
correlated in x and H. The LIM corresponding to the Manna model has a noise field rj(x, H) 
which is point-like and correlated: in the iJ-direction because of random- walk like increments 
in r](x, H) and in the x-direction because of the short-ranged grain conservation (if a nn of x 
gained a grain when H(x, t) increased other nn's did not) 

For non-SOC (periodic or open) boundary conditions the LIM has a depinning transition 
at a critical force F c . At the transition, the scaling exponents of the LIM with columnar, 
point-like, and, depending on the details, correlated disorder differ (see ||,H] for the effect of 
noise correlations). This means that avalanches have different spatial and temporal properties 
since the critical exponents like the roughness exponent \ of the LIM depend on the noise. 
In particular, the BTW model is in a different universality class from the others as 
it has no 77-noise. The LIM is invariant to forces that are static in the 7J-direction ||| which 
makes r\ a relevant perturbation. For the Manna-model, we conjecture that it may be in the 
point-disorder LIM class despite the correlations in 77 0,0). This prediction seems to be shown 
to be true in d = 2 in Ref. p5| . The upper critical dimension of the LIM is d u = 4 for all 
these kinds of noise because of the Laplacian in Eq. (||) . 

These conclusions are based on the continuum limit of the LIM. However, sandpile models 
are discretized versions thereof, with the additional cr-noise term in the equation. In the point- 
disorder LIM, various numerical approaches indicate that the a is irrelevant for the rice-pile 
model with periodic boundary conditions This can be understood by considering the 

sum of the 77- and cr-noise terms: the velocity constraint just renormalizes the 77-field since 
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the velocity would in any case be of the order of unity, or, the actual value of cr depends on 
z c (x,H). 

For the other scenarios one has to answer the question, whether the cr-noise changes the 
universality class of the model at hand. This would mean simply that the sandpile model can 
be a "bad discretization" of a continuum interface equation so that the cr-noise changes the 
properties of the continuum model. Notice that the way o~(x, H) becomes non-zero (in these 
models a < 0) depends on the model. In the Manna-model a site with z x = 1 can also get 
two grains from the same nn leading to a ^ 0. This is point-like disorder since the decision 
to give two grains has no correlations with any earlier events. Thus the Manna model may 
be in the rice-pile universality class as noted above. Like in the rice-pile model, the cr-noise 
should not have any strong correlations due to the randomness in the avalanches. 

In the BTW and rice-pile models the cr-term arises (see again Fig. |l|)) from the simultaneous 
topplings of nn's of site j. Thus in d > 2 the role of the noise is diminished, and in particular it 
should not affect the upper critical dimension. Notice that the cr-noise changes off the critical 
point (since more of the neighbors are likely to be active at the same time) , and may thus play 
a role in systems off criticality like in the 'fixed density' (or energy) ensemble corresponding 
to a constant total force fTo| . The BTW model has also the Abelian property || and thus the 
cr-field of any particular sample depends on the exact order of the topplings. 

The cr-noise is next studied as such to elucidate the role of the SOC boundary conditions 
in the interface equation (open boundaries for the sandpile). We look at (mostly) the BTW 
and the other models by numerical simulations in 2d in the normal SOC sandpile ensemble 
using parallel dynamics for all active sites at each time step, as one would do in an interface 
model. Quantities of interest are the average of cr, and its spatial dependence on x. The 
cr-field is constructed from the relation (|]). 

First we look at the finite size dependence of the probability to have a non-zero a, Pl(o~ < 
0). It can be studied as averaged over dissipating (grains leave the system), and non-dissipating 
avalanches. For all these Pl increases with the system size L. The finite size scaling Ansatz 
Pl ~ (-P)oo + aL~ b , fitted to the data |l6| ], gives correction exponents that are close to 
6 = 2/3 for bulk and dissipating avalanches, and —1 for all avalanches; a is negative in all 
cases. The asymptotic BTW values are (P) ~ 0.081 for all and (P) d — 0.121 for dissipating 
avalanches. Thus the dissipating avalanches have typically stronger cr-noise: the interface tries 
to move faster. They also contribute significantly into (cr) since the total number of topplings 
is dominated by such avalanches. Therefore fluctuations of the cr-field may be related to the 
fluctuations in the average grain number, (z x ) ~ { z x) c - For the Manna and ricepile models 
(a) also decreases with L [Q. 

We checked the spatial correlations in o~(x, H) by computing the two-point noise-noise cor- 
relation function C(5H) = (a(x, H + SH)a(x, H)} — (cr 2 ), at a fixed x — const. [0. Figure || 
demonstrates that at the center and at the edges of the pile the noise decays exponentially. 
The decay is slower in the center the decay lengths being in general proportional to L. In the 
'bulk' (see the second and third curve from left) we see superposed on that behavior periodic 
oscillations. Thus there is no spatial translational invariance |p"8|-po|| since the correlations 
depend on x. The details of the correlations of cr in the BTW model may be related to ob- 
servations of multiscaling in the BTW model pl| and manifest the columnar character of the 
F-noise. In contrast, for the Manna/rice-pile models the correlations in cr decay rapidly. 

The average (o~(x)) turns out to be non-uniform in x. For instance, at the boundary 
the noise is weaker since sites can receive grains only from 2d — 1 nn's. The non-uniform 
cr implies simply that the interface tries to move faster in the center of the pile. Next we 
investigate whether a(x, H) ^ if site x does not topple any more during an avalanche. In 
other words, is the velocity constraint related to the stopping properties of avalanches? We 
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look thus at the cr-noise at the elastic pinning paths of an interface model ]22| . The answer is 
relevant for the existence of translational invariance as well. Figure || shows a data collapse 
of P(a(x,y) ^ 0) along the cut y = L/2, 1 < x < L/2, scaled with P(x ~ L/2), with the 
pinning-path constraint, i.e. the fraction of all topplings with a < such that the toppling 
is the last at the site during an avalanche. The probability increases at the boundaries with 
system size L. This shows that the standard BTW sandpile has no spatial translational 
invariance. This is true (see the inset) also for the Manna and rice-pile models. Their pinning 
paths are determined by the configuration prior to the avalanche and the (point) disorder 
the avalanche encounters. In the BTW case the paths are set by the original force/grain 
configuration as in a columnar-noise LIM, but not by the cr-noise field itself as the Abelian 
character of the BTW also implies (the choice of the toppling order does not change the final 
configuration) . It would be interesting to compare the BTW result to the wave picture of the 
BTW model |^|,|2^,|2j]. The lack of translational invariance is caused by the open boundary 
conditions that also imply a parabolic shape for the interface and suggests that there is no 
simple scaling as one would expect for a normal interface model, in which case a relation 
l D ~ l d+x would be valid for the avalanche size vs. its linear dimension I. At most, one should 
have D s < d + x for the cut-off dimension of the probability distribution of the avalanche 
sizes, since the open sandpile avalanches can of course not be 'over-critical' with any effective 
x' > X- Note that this observation seems to be true for any of the three models. 

In conclusion, sandpiles can be mapped to driven interfaces, by describing the dynamics 
with various types of quenched noise. The rice-pile model is equivalent to the random field 
linear interface model. The BTW model has long-range on-site correlations that can be 
studied via the cr-noise or the restriction v < 1. It is because of the deposition noise F(x, t) a 
columnar-noise LIM albeit with the velocity limitation. The Manna model turns out to have 
"correlated point-disorder" and is thus in a different universality class from the BTW model. 
The projection technique used for the Manna toppling dynamics can be applied to e.g. the 
Zhang model, models with bulk dissipation of grains, and the Olami-Feder-Christensen model 
p5|-p7t and thus further extensions of our work are certainly possible. The discussion of the 
cr-term in the interface equation for sandpiles makes it clear that studying the noise provides 
a new tool for elucidating sandpile behavior and the role of various boundary conditions in 
sandpiles. The fundamental ingredients in sandpiles are slow drive and fast dissipation. With 
open boundaries these combine to make the pile spatially non-uniform, while the rules chosen 
are reflected in the scaling of the self-organized critical state. 

* * * 
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Fig. 1 - Rescaling of the force / and an example of how the a-noise ensues (three grains added 
simultaneously) . 
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Fig. 2 Fig. 3 

Fig. 2 - Plots of C(Ah) for BTW model noise at locations ranging from the edge to the center (left 
to right) for L = 64. Inset: C(Ah) in the center for L = 64, 128, 256, 512 (bottom to top). 

Fig. 3 - P(cr(x) < 0) vs. x scaled with P(x = L/2) for various L in the BTW model (L = 32 (circles), 
64 (triangles), 128 (squares). Inset: for L = 64 the same for BTW, Manna, ricepile (triangles, circles, 
crosses) . 



